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Introduction. Ising quantum chains with aperiodi- 
cally modulated coupling constants have been the topic 
of several recent publications,!!! - !] see also papers cited in 
Ref. 0. These systems describe 2d classical Ising models, 
with couplings varying along layers, in the Hamiltonian 
limit and provide a non-trivial, but nevertheless easily 
accessible class of models to study the effect of aperiodic 
(dis)order on thermodynamic properties and critical be- 
havior. Scaling arguments, put forward by Luck, lead 
to the following predictions on the relevance of aperiodic 
modulations □ The critical behavior is of Onsager-type as 
long as the fluctuations of the couplings are bounded, but 
resembles the randomly disordered case for unbounded 
fluctuations. Of particular interest is the marginal case 
with (in terms of the system size) logarithmically diverg- 
ing fluctuations. Here, a non-universal scaling behavior 
is predicted with critical exponents that depend contin- 
uously on the coupling constants. 

For sequences of coupling constants obtained from sub- 
stitution rules, we demonstrated recentlyjj how Luck's 
criterion can be derived within an exact renormalization 
scheme, using a decimation procedure proposed in Ref. ^. 
This way, an exact expression for the scaling exponent z 
of the mass gap at criticality was obtained. 

In this article, we show how the renormalization ap- 
proach can be extended to describe the fermionic eigen- 
vectors of the quantum chains. This leads to the exact 
determination of the critical scaling behavior of several 
bulk and surface quantities, like the bulk energy density 
or the surface magnetization. Quite recently, a similar 
derivation (using the decimation of Ref. [|) has been given 
for some special examples.B These results are contained 
as special cases here. 

Ising quantum chains and substitution sequences. 
The Ising quantum chain is given by the Hamiltonian 



H 



N 




(1) 



with Pauli matrices cr x k ,z 



acting on the fcth site, and 
site-dependent transverse magnetic fields hk- The site- 
dependent coupling constants Ek shall be drawn from a 



set of values e ai , where the label ai is taken from an n- 
letter alphabet A = {a±, 0,2, . . . , a n }. Initially, the field 
variables are chosen as hk = 1; but later on, as param- 
eters in the renormalization transformation (RT), they 
become site-dependent and will correspondingly be la- 
beled as hk+i = h ai (if e k = Sen)- 

We choose the couplings according to sequences gen- 
erated by iterated application of a substitution rule 
q : ai — » w(ai) on A. Some important properties of the 
substitution chain are already contained in the corre- 
sponding n x n substitution matrix M B , whose elements 
(Mg)^ give the number of letters aj contained in the 
words w(a,i). So, the asymptotic densities of the letters 
are given by the entries of the (statistically normalized) 
right eigenvector to the Perron-Frobenius (PF) eigen- 
value Ai = App which is just the rescaling factor of the 
chain length. The fluctuations of the cumulated devia- 
tion of the couplings from the mean coupling grow as a 
power of the chain length, i.e. as N u , with a wandering 
exponent u) determined by the two largest eigenvalues (in 
modulus) 



u> = 



log |A 2 | 
log Apf 



(2) 



Thus fluctuations diverge for |Aa| > 1, whereas the 
marginal case is connected to |A 2 | = 1. 

The Ising quantum chain is critical, in the sense that 
the mass gap vanishes in the thermodynamic limit, if 
the geometric mean of the reduced couplings Sj/hj ap- 
proaches 1 for N — > 00. For couplings determined by 
substitution sequences, this is equivalent to the condition 
that the vector y of logarithms of the reduced couplings 



y, 



log 



(* = 1,2, 



(3) 



is perpendicular to the PF eigenvector of M e . 

For a general set of coupling constants e k , the quan- 
tum chain can be written as a free-fermion model via a 
Jordan- Wigner transformation and diagonalized canoni- 
cally, resulting inEI 
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where rj q and rj q are N fermionic creation and annihi- 
lation operators and where CV is some constant. The 
dimensionless excitation energies A ? (which can be or- 
dered as < Ai < A2 < . . . < Ajy) satisfy the linear 
difference equations 



(5a) 
(5b) 



This corresponds to an eigenvalue problem of the matrix 
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where e N = for free boundary conditions which are 

employed in the study of surface properties. Here, ip^ 

and 0[ are the even and the odd components of the 
eigenvector connected to the qth eigenvalue, respectively. 

Renormalization scheme for fermionic eigenvectors. 
Following the presentation in Rcf. ^, we now introduce 
the renormalization technique,based on <S-transfer matri- 
ces and their star-products.E3 For technical reasons, we 
restrict the discussion to substitution rules of the form 



(7) 



where Wi are n arbitrary words in the n-letter alphabet 
A, and only comment on the extension to the general 
case. We introduce 5-matrices that transform the com- 
ponents of the eigenvectors of Eq. (||) as 



k\I 



( 4>k+i 
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with 



$k\l = <5 fe | fe+1 * S k+l \ k+2 * . ..* Si_i\i , (9) 
where the *-product of two 2x2 matrices is defined as 
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According to Eq. (]6[), the elementary «S-matrices are 



'k\k+l 
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where A denotes an eigenvalue of 7i. In order to estab- 
lish the RT, the fields and vector components arc labeled 
according to the locally varying coupling constants. Fur- 
thermore, two additional asymmetry coefficients k^. have 
to be introduced for every letter ai € A. This way, we ob- 
tain n 2 different elementary <S-transfer matrices through 



^ai,k \ _ 



4>ai,k+l 



(12) 



^4>a,j ,k+2 , 

where the second label refers to the position, and where 



& ai I a j 



(13) 



Unlike the fields or the couplings, the eigenvector com- 
ponents depend not only through the a, , but also explic- 
itly on the position, so we have to keep the labels here. 
The RT now reverses the substitution steps by building 
^-products of <S-transfer matrices corresponding to the 
words Wi given by the substitution rule. The resulting 
RT's of kJ. A and h ai have been described in Ref. M. For 
coupling constants fulfilling the criticality condition, the 
vector y of logarithms of the reduced couplings is found 
to scale with the next-to- leading eigenvalue A2 of M e . 
Especially, in the marginal case with A2 = 1, y converges 
to an eigenvector corresponding to an eigenvalue 1 of the 
transpose M l g of the substitution matrix. 

Let us now concentrate on the RT equations for the 
vector components of ij) and <p. Since iterated *-multi- 
plications of the «S-matrices act as a pure decimation on 
the vector components, the remaining components of the 
renormalized vectors tp and <f> differ from the original 
ones only through relabeling and a normalization factor 



tpm = C^lpk , 0m+l =C^4> 



fc+1 



(14) 



with k/m — > Apf and normalization constants ^ such 
that both if) and <f> are normalized to one. We now pro- 
ceed to determine the scaling behavior of these normal- 
ization constants at the critical point A = 0. As the 
diagonal elements of the <S-matrices vanish to leading or- 
der in A, the RT equations for tp and <fi decouple at the 
critical point and can be considered separately. By tak- 
ing successive ^-products along the words Wi , the original 
eigenvector components may now be expressed in terms 
of the renormalized ones 



Tpa t ,k 



'Pa^k+l 



where w\ denotes the jth letter of Wi and 

S a z \wl = S a t \w\ * S wl\w* S wj- 1 \ 

Explicitly, we obtain the following relations 
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with the reduced couplings at their fixed point values. 
In order to derive the renormalization factor of C$ and 
C^, we now perform partial sums of the squared vector 
entries, according to their additional labels a,, 



(19) 



For these vectors ^ and the RT at the critical point 
yields a simple matrix form 



C|* = M v '* , Cl$ = M I> , 



where 
M 



(20) 
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Since all components of the vectors $ and the ma- 
trices are positive, the vectors converge to the PF 
eigenvectors of M^>^ under iteration of the RT. We for- 
mally conclude, for the normalization constants, 



G,i, — 



(22) 



where are the PF eigenvalues of M^^. At this 

point, we can make contact to the RT of the low- 
est fermionic excitations considered in Ref. [| To this 
end, we transform M^'^ under conservation of the spec- 
trum according to M = (T^)" 1 M^T^ and M+ = 
(T ) -1 M^T^ with diagonal transformation matrices 



T% = (e a Jh ai ) 2 S hJ and T% = h 2 <),.,. We find 



|ti>,| k / h 2 



±1 



(23) 



k=l 
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If /i = fi^ and /i + = yfi denote the leading eigenvalues 
of M ± , note that /i~ '({e a Jh ai }) = ^ + ({/i a ,/£a,})- 

For substitution chains with bounded or marginal fluc- 
tuations of the coupling constants, however, these are 
exactly the forms of the RT-matrices for the fermion fre- 
quencies [see Eqs. (3.21) and (3.25) in Ref. |[. For com- 
parison, recall that one has 
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(24) 



for | Aa j < I- The lowest fermion frequencies transform 
as A = yj /i+/i~A, resulting in the scaling behaviorH 

log(M + M") 



A„ 



x q N~ z 



2 log Apf 



(25) 



The extension of the RT to general substitution rules 
has been described in Ref. |[ it can be performed along 
the same lines here. Generally, the substitution rule has 
to be redefined as a function acting on pairs of letters, 
generating the same chain. General «S-matrices S CH \ aj 
contain fields and asymmetry parameters, and act on vec- 
tor components, all carrying a double label (. . .) ot aj ■ The 
RT of the eigenvector components can again be given as 
a matrix equation (of dimension n 2 ) and identified with 
the RT of the fermionic excitations as above. 

Results on critical scaling. The surface magnetization 
m s i (on the left surface) can be obtained from the large 
distance limit of the spin-correlation function at the sur- 
face, leading to the expressionErtZl 



(lK|O)=0i 



(i) 



(26) 



Here, |1) = rf\\0) is the lowest excited state which be- 
comeS|-degenerate with the ground state in the ordered 
phased and the last equality is obtained upon express- 
ing af in terms of the fermions. Thus, m s i is simply 
the first entry of the normalized eigenvector (jy- 1 ' of the 



lowest fermionic excitation. Since <p^ = yj fi + (j)\ 



<>(!) 



+ = 



t 1 



[see 
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Eqs. ( |l4| ) and (p2|), and remember that \x 
obtain the surface magnetization exponent (3 s i through 



m s i 



AT 
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2 log Apf 
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Note that [3 s i is given by the scaling exponent here, since 
the correlation length exponent (along the layers of the,2d 
classical model) retains its unperturbed value v = 1J10 
Using a method of Ref. |lj, the surface magnetizatioa 
of several substitution sequences was calculated before.EJ 
However, it seems impossible to extend this method to 
the general casg. Recently, the results of Ref. [j] have also 
been rederivedfj using real space renormalization. 

For substitution chains with bounded fluctuations (and 
second largest eigenvalue of the substitution matrix 
| As | < 1), like e.g. the Fibonacci chain, all reduced cou- 
plings converge to the critical value of the uniform chain 
ZaJKi -> 1 and we find M+ = M = M e , thus 
H + = Xpf and [3 S \ takes its Onsager value (3 S \ = 1/2. 

For two-letter substitution chains, which alwayS|-|Can 
be generated by a substitution rule of the form (fn) p an 
explicit result can be obtained also for the marginal case 
(A2 = 1). Here, the PF eigenvalue of coincides with 
the trace of these matricesjj resulting in 



log(M+ + M+-1) 



2 log A 



PF 



(28) 



where Mfi and M%2 are given in Eq. (|2^). In general, 
[3 si depends on the coupling constants through r < n 



3 



independent parameters, where r is the dimension of the 
(joined) eigenspace to eigenvalues |A| = 1 of M e . 

For substitution chains with fluctuations that diverge 
with a power-law, the situation is slightly more com- 
plex. Here, the reduced couplings scale as exp^jiV^); 
and since y always contains components of different sign, 
part of them flow to infinity under the RT, while oth- 
ers converge to zero. For the RT of C^, consider the 
behavior of the entries of 7W ± for iterated substitution 
steps g n . According to Eq. (|23|), all entries are polyno- 
mials of the reduced couplings. Now, two cases have to 
be distinguished — depending not only on the substitu- 
tion matrix, but on the detailed form of the substitution 
itself. Generically, all entries are finally dominated by re- 
duced couplings that flow to infinity under the RT. This 
is always the case if the cumulated deviation of any type 
of coupling from its mean is unbounded in positive and 
negative direction. In this case, m s \ will finally scale as 
exp(— cN^), with c depending on the coupling constants. 

There are, however, special cases, where the cumulated 
deviation does not diverge symmetrically, but, depending 
on the initial word (seed) of the inflation, only in positive 
(or in negative) direction. As a simple example, consider 
the two-letter substitution rule 



Q ■ 



bbba 



(29) 



Let §a,b(n) be the number of a, 6's among the first n let- 
ters of the limit chain; we find 1/2 as the asymptotic 
frequency for both letters. Choosing, however, the sin- 
gle letter a as the seed, the number tt a (n) — %(n) is un- 
bounded, but positive for any n, thus the fluctuations are 
bounded from below for this chain. For the correspond- 
ing matrices M^, this means that, for given couplings, 
certain entries do not diverge, but finally tend to a finite 
value or to zero. Starting with the appropriate seed, this 
results in a finite surface magnetization m s \ (on one side 
of the chain), also at criticality. This explains the (un- 
usual) behavior that has been observed in Ref. || for a 
special substitution rule. 

The surface magnetization m S 2 on the right surface 
follows, after an inversion of the substitution rule to 
g : ai — > WiOi, through the same procedure as described 
above. The RT for the vector components reads 



' = s 



w i \ a i 



(30) 



and we deduce the same relations as above with all re- 
duced couplings inverted, i.e., e a Jh ai — -> h a Je ai - Thus 

^({? 1 })=^({r i }) = o^- (3D 

fl ai £ a% £ log Apf 

This generalizes a symmetry in scaling of m s \ and m^&^l 
also to substitution chains where this symmetry is not 
obvious from the chain geometry, and to general n-letter 
substitutions. We also conclude, for arbitrary n-letter 



substitution chains with marginal (or bounded) fluctua- 
tions of the coupling constants, that 

z = 0„i+ /3 s2 (32) 

as conjectured in Ref. || and argued-for in Ref. [g. 

The surface energy density can be obtained from the 
two lowest fermionic excitations as followsa 

Ai) 



e«i 



(A 2 



(33) 



Hence, we find the scaling behavior e s \ : 2 ~ N~ z ~ 2 P s1 - 2 
conjectured in Ref. ^ on the basis of numerical results. 

The bulk energy density on site j follows from the ma- 
trix elemental 



eft 



v>iM 2) -v> fc 



N~ 



(34) 



as had also 



Thus, e m = y/n+n e k and e N/k 
been conjectured in Ref. ||[ 

Conclusions. We have demonstrated how the bulk 
and surface critical behavior of aperiodic Ising quantum 
chains (or, equivalently, of aperiodically layered, two- 
dimensional classical Ising systems) obtained from gen- 
eral substitution rules can be derived by exact real space 
renormalization. The results corroborate a number of 
conjectures on scaling laws for these systems and gener- 
alize previous results obtained for special examples. An 
exact calculation of the bulk magnetization seems to be 
more difficult and remains an open problem. 
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